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ABSTRACT 

An optically and geometrically thick torus obscures the central engine of Ac- 
tive Galactic Nuclei (AGN) from some lines of sight. From a magnetohydrodynamical 
framework, the torus can be considered to be a particular region of clouds surrounding 
the central engine where the clouds are dusty and optically thick. In this framework, 
the magnetic field plays an important role in the creation, morphology and evolution 
of the torus. If the dust grains within the clouds are assumed to be aligned by para- 
magnetic alignment, then the ratio of the intrinsic polarisation and visual extinction, 
P(%)/Av, is a function of the magnetic field strength. 

To estimate the visual extinction through the torus and constrain the polarisation 
mechanisms in the nucleus of AGN, we developed a polarisation model to fit both the 
total and polarised flux in a 1.2" (~ 263 pc) aperture of the type 2 AGN, IC5063. The 
polarisation model is consistent with the nuclear polarisation observed at Kn (2.0 - 
2.3 |xm) being produced by dichroic absorption from aligned dust grains with a visual 
extinction through the torus of 48±2 mag. We estimated the intrinsic polarisation 
arising from dichroic absorption to be P^"^ = 12.5±2.7%. 

We consider the physical conditions and environment of the gas and dust for the 
torus of IC5063. Then, through paramagnetic alignment, we estimate a magnetic field 
strength in the range of 12 - 128 mG in the NIR emitting regions of the torus of 
IC5063. Alternatively, we estimate the magnetic field strength in the plane of the 
sky using the Chandrasekhar-Fermi method. The minimum magnetic field strength 
in the plane of the sky is estimated to be 13 and 41 mG depending of the conditions 
within the torus of IC5063. These techniques afford the chance to make a survey of 
AGN, to investigate the effects of magnetic field strength on the torus, accretion, and 
interaction to the host galaxy. 

Key words: AGN, torus - infrared: polarimetry. 
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The detection of polarised broad emission lines in the spec- 
trum of NGC1068 revealed an obscured Active Galactic Nu- 
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cleus (AGN) through scattering of the broad line region 
(BLR) radiation (Antonucci fc Miller|1985[|. This study gave 
a major boost to the unified model ( Antonucci||1993 Urry 



fc Padovani|1995| ) of AGN, which hods that all AGN are es- 
sentially the same object, viewed from different line of sight 
(LOS). In the unified model scheme, the AGN classification 
solely depends on the obscuration of an optically and geo- 
metrically thick dusty torus 



Recent studies (e.g. Nenkova et al. 2002 2008 1 have 



proposed that the dust in the torus is distributed in clumps, 
and accounts for (a) the distribution of temperatures along 
the torus; and (2) the variety of spectral energy distribu- 
tions (SED) produced by geometry, clumpy distribution and 
spectral features, such as the 10 \im absorption and emis- 
sion. The clumpy torus model holds that the dusty clumps 
are distributed in a few parsecs, consistent with observa- 
tions (e.g. Radomski et al.|2002 Jaffe et al.|2004 Packham 



et al.|2005[|Suganuma et al.|2006| [Tristram et al.|2007[ ). This 
means that the torus cannot be resolved at optical/infared 
(IR) wavelengths, even with the high-spatial resolution pro- 
vided by 8-m class telescopes. Recent high-spatial resolution 
observations of Seyfert galaxies have shown that the clumpy 
torus model can account for the near-IR (NIR) and mid- 
IR (MIR) emission (see 'Mason ot al. 2006 , Ramos Almeida 
et al.||2009[ |2011| [Alonso-Herrero et al.,,2011). Some stud- 
ies ( |Mor, Netzer, fc Elitzur|2009| |Mor fc Netzer|2012| ) have 
shown that a component composed by hot graphite dust can 
explain the NIR emission of Type 1 AGN. The clumpy torus 
models permit, from a statistical view, an examination of the 
general properties of the torus, i.e. inclination to our LOS, 
number of clouds, covering factor, optical depth of clouds, 
etc. However, the intrinsic properties of the clumpy torus, 
i.e. the dust grain composition, grain size, grain alignment, 
etc. remain unknown. Polarimetry techniques show a pow- 
erful potential to obtain the intrinsic properties of the torus. 



Several polarimetric studies in the IR to NGC1068 (Mason 



et al. 2007 Packham et al. 20071 have demonstrated this 
potential, constraining (1) interstellar dust properties; (2) 
upper-limit diameter of the torus; (3) mechanisms of polari- 
sation and hence magnetic field directions in the central few 
pc. 

Several models to explain the existence of the torus 
have been proposed. Some models explain the torus as an 
infiow of gas from large scales (e.g. Wada et al.|2009 Schart- 



mann et al.||2011 l. Wada et al. (20091 presented numerical 



simulations of the interstellar medium to track the forma- 
tion of molecular hydrogen forming an inhomogeneous disc 
around the central engine, identified as the torus. |Schart-| 
|mann, Krause, fc Burkert] ( |2011[ ) suggested the origin of the 
torus as clouds falling to the central engine. In this scheme, 
it is difficult to explain how the energy dissipation of clouds 
and cloud-cloud collision tends to concentrate the clouds 



to form the torus near the midplane (Krolik & Begelman 



19881. Conversel; 


f, some models (Blandford & Payne||1982 


Emmering et al. 


1992 Elitzur fc Shlosman 20061 suggest 



that the clouds are confined by the magnetic field generated 
in the central engine and are accelerated by the hydromag- 
netic wind. In this scheme, the hydromagnetic wind could lift 
the clouds from the midplane to form a geometrically thick 
distribution of clouds surrounding the central engine, where 
the torus is in a particular region of the hydromagnetic wind 
where the clouds are dusty and optically thick. |Kartje et al.| 



( 1999 1 suggests that magnetic field strengths greater than 20 
mG can account for a clumpy disc-driven wind model, where 
clouds moves along the magnetic field lines in a homogeneous 
outflow component. VLA circular polarimetry observations 
of NGC4258 inferred an upper-limit of the magnetic field 
strengths of 300 mG 



fHerrnstcin ct al. 1998). Further po- 
larimetric observations (Modjaz ct al. 2005) at 22GHz of 
the water vapour maser clouds in NGC4258 using the VLA 
and GBT, estimate magnetic field strength from 90 to 300 
mG at a radius of 0.2 pc from the central engine. In sum- 
mary, observations and theoretical models for maser clouds 
at distances ~0.2 pc seems to be in agreement. However, ob- 
servational or theoretical studies at larger distances, where 
the clouds are optically thick and dusty, have not been done 

to date. 

The nearby (z = 0.0]p| |de Vaucouleurs et"al|(|l99l|) e l- 



liptical galaxy (r^^** brightness profile; Colina et al. (1991l) 
IC5063 (PKS 2048-572) shows polarised scattered He broad 



emission lines suggesting a hidden Type 1 AGN ( Inglis et al. 
|1993j ), and an increase in the degree of polarisation in the 
infrared at J, H and K (Hough et al. 19871. A prominent 



dust lane has been observed ( Colina, Sparks, & Macchetto 
|1991[ ) along the long-axis of IC5063. IC5063 shows a ra- 
dio lumino sity (log Piaghz = 23.8 W Hz-\ Ho = 50 km 



s kpc , Colina, Sparks, fc Macchetto (1991l) about two 



orders of magnitude larger than is typical nearby Seyfert 
galaxies, classify IC5063 in the range of low-luminosity ra- 
dio elliptical galaxies. The central polarised source have been 
suggested to be a BL Lac object from infrared observations 



by Hough et al. (19871. They showed a steep-spectrum IR 
component and suggested a synchrotron central source. Fur- 
ther NIR studies at J, H and K bands by |Brindle et al.| 



(|1990 1 modeled the total and polarised flux of IC5063, and 
concluded that the nuclear source can be explained by a red- 
dened central source with an power-law index of 1.5 (typical 
for Seyfert 1). 

In this paper we estimate the magnetic fleld strength 
in the NIR emitting regions of the torus of AGN, through 
NIR polarimetric observations of IC5063 in the J, H and 
Kn bands. We develop a polarisation model to simultane- 
ously fit the total and polarised flux, that allow us to inter- 
pret the different mechanisms of polarisation to the central 
engine of IC5063. Further estimates of the extinction to the 
central engine at different wavelengths allow us to interpret 
the origin of the NIR polarisation. Finally, four independent 
methods of calculating the magnetic field strength in the 
clumpy torus of IC5063 are presented. 

The paper is structured as follow. In Section 2, we de- 
scribe the observations and data reduction. In Section 3, 
the results are shown, then analysed in Section 4. We dis- 
cuss various mechanisms of polarisation that could account 
for the polarisation in the nucleus of IC5063 and we develop 
a polarisation model to interpret our data in Section 5. In 
Section 6, the nuclear extinction estimated at various wave- 
lengths is discussed. In Section 7, the magnetic fleld strength 
of the torus of IC5063 is estimated. Finally, conclusions are 
presented in Section 8. 



^ Through this work, we adopt Hq = 73 km s ^ kpc ^ , so that 
1" = 219 pc for the redshift of IC5063 



© 2012 RAS, MNRAS 000,[T|fT6l 



Magnetic field strength in the torus of AGN 3 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made using the inharcd polarimeter 
built by the University of Hertfordshire (Hough ct al.|1994 1 
for the Anglo- Austrahan Telescope (AAT). The polarimeter 
features consist of an achromatic (1 - 2.5 yuai) half- wave 
plate (HWP) retarder from 1.0 to 2.5 \im, which can be 
stepped to set angular positions. This is placed upstream of 
the observatory's near-infrared camera, IRIS ( Gillingham &| 



Table 1. Summary of observations. 



Lankshear| 1990 1 , in which a magnesium fluoride WoUaston 
prism was used as an analyser. A mask in the focal plane of 
IRIS served to blank out half of the field so that the ordinary 
(o-ray) and extraordinary (e-ray) ray from WoUaston prism 
do not overlap when imaging extended objects as well as 
reduce the sky from overlapping. 

IC5063 was observed on 1994 July 25, in the J, H and 
Kn (2.0 - 2.3 Jim) bands, with IRIS at the f/15 focus using 
the Rockwell camera with 128 x 128 pixels. The pixel-scale 
is 0.6" per pixel, giving a field of view of 76" x 76". The 
half-wave retarder was stepped sequentially to four position 
angles (PA; 0°, 45°, 22?5 and 67?5) taking an exposure at 
each HWP PA. Images at each HWP PA were flat-fielded, 
sky-subtracted and then cleaned by interpolating over dead 
and hot pixels and cosmic rays. Next, the images were reg- 
istered and shifted by fractional pixel amounts in order to 
account for slight image drift between frames. Then, the po- 
larisation images were constructed using our own IDL rou- 
tines. The individual o- and e-ray images were co-averaged 
to increase the signal-to-noise ratio. Next, the o- and e-rays 
were extracted using a rectangular aperture of 28" x 18". 
The Stokes parameters, I, Q and U, were calculated accord- 



ing to the ratio method prescription (see Tinbergen ( 19961) 



Finally, the degree of polarisation, P, was derived such as: 



P = 100 X 



v/Q2 + 1/" 



and the position angle of polarisation, 6, was found by: 



tan 



A summary of observations is shown in Table [T] The 
night was photometric with the seeing estimated to be 1.2 
from the automatic-guide at the observatory. Flux standard 
stars were not observed during this night, so we used the 
photometric data from Axon et al. ( 1982 I . Specifically, the 



flux calibration was performed using photometry at J, H and 
K in a 10" aperture, where the large aperture ensures that 
centering issues are minimal. The difference between the K 
and Kn bands makes a small contribution to the total photo- 



metric error. The zero-flux standard values ( Campins et al 



1985[ ) of 1603 Jy, 1075 Jy and 667 Jy for J, H and K„, re- 
spectively were used to transform between magnitudes and 
fluxes. To calibrate the PA of polarisation at J, H and Kn, a 
4.5" aperture from Hough et al. (19871 was used. Since un- 
polarised standard stars were unavailable, an instrumental 
polarisation of 0.02±0.06%, from previous studies by |Hough| 
et al. ( 1994 \ was adopted. 



Filter 


Frame Time Sets 


Total observation time 










(s) 


J 


30 


12 




1440 


H 


20 


8 




640 


Kn 


10 


8 




320 


Table 3. The observed degree and position 




tion and photometry of IC5063. 






Aperture Filter 


P 


e 


Magnitude 


(") 




(%) 


n 




1.2 


J 


2.0 ± 0.7 


± 2 


14.81 ± 0.02 




H 


2.5 ± 0.9 


5 ± 8 


13.90 ± 0.02 




Kn 


7.8 ± 0.5 


4 ± 4 


13.01 ± 0.06 


2.0 


J 


1.5 ± 0.5 


-1 ± 2 


13.88 ± 0.02 




H 


2.0 ± 1.6 


4 ± 11 


12.95 ± 0.02 




Kn 


6.2 ± 0.3 


4 ± 3 


12.26 ± 0.03 


3.0 


J 


1.1 ± 0.3 


2 ± 4 


13.28 ± 0.01 




H 


1.6 ± 0.9 


3 ± 7 


12.38 ± 0.01 




Kn 


4.7 ± 0.3 


4 ± 2 


11.79 ± 0.02 


4.0 


J 


0.9 ± 0.2 


5 ± 5 


12.89 ± 0.01 




H 


1.1 ± 0.5 


3 ± 6 


12.01 ± 0.01 




Kn 


3.8 ± 0.1 


4 ± 2 


11.48 ± 0.01 



3 RESULTS 
3.1 Photometry 

We made measurements of the nuclear total flux in several 
apertures to compare to previously published values (Ta- 
ble |2|. In aU cases, photometric errors were estimated by 
the variation of the counts in subsets of the data. The dif- 
ference in fluxes between our and previously published re- 
sults is < 4% in all apertures. To optimally measure the 
AGN flux of IC5063 and minimize the contributions from 
the host galaxy, dust lane and ionisation cones, photometry 
in an aperture equal to the seeing of the observations (1.2") 
was used. To investigate the flux dependence on the aper- 
ture size, photometry in three apertures, 2.0", 3.0" and 4.0", 
was measured (Table [3|. An increase in total flux at both 
longer wavelengths and larger aperture was found (Figure 
[l]). In line with previous studies, an elliptical profile with a 
FWHM of 3.7" X 2.6", 3.4" x 2.5" and 2.3" x 1.8" at J, 
H and Kn respectively, was found for the nucleus of IC5063. 
In all filters, IC5063 was found to be extended along PA — 
75°. The FWHM was estimated using a Gaussian profile. 



3.2 Polarimetry 

We made measurements of the nuclear degree of polarisa- 
tion in several apertures to compare to previously published 
values (Table [2|. In aU cases, polarimetric errors were esti- 
mated by the variation of the counts in subsets of the data 
The observed degree and PA of polarisation in four different 
sized apertures, 1.2", 2.0", 3.0", 4.0", were measured (Ta- 
ble |3| . Polarimetry maps of the total and polarised flux are 
shown in Figures [2] and [3] in J, H and Kn bands. In these 
figures, the overlaid polarisation vectors are proportional in 
length to the degree of polarisation with their orientation 
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Table 2. Comparison with literature of the nuclear photometry of IC5063. 



Aperture Filter 
(") 



Hough et al. (19871 
Mag P 
(mag) (%) 



This work 
Mag P 
(mag) (%) 



2.25 
4.5 



H 

Kn 
J 

H 

Kn 



12.38 
11.79 
12.50 
11.80 
11.10 



1.75±0.34 
6.33±0.50 
0.77±0.07 
0.99±0.05 
3.25±0.45 



12.86 ± 0.02 
12.18 ± 0.02 
12.71 ± 0.01 
11.83 ± 0.01 
11.34 ± 0.01 



1.8±1.5 
6.0±0.3 
0.7±0.2 
1.0±0.4 
3.3±0.1 



Aperture Filter 
(") 



Axon et al. ( 1982|) 
KTag 
(mag) 



This work 
Mag 
(mag) 



5.0 



J 
H 



12.72 
11.80 
11.32 



12.59 ± 0.01 
11.72 ± 0.01 
11.23 ± 0.01 




Wavelength (/im) 



) 2.0 

elength (/im) 





Wavelength (^im) 



1,5 2.0 
Wavelength (^m) 



Figure 1. Total flux {top left), degree {top right) and position 
angle {bottom left) of polarisation and polarised flux (bottom 
right) versus wavelength. Measurements with aperture size of 1.2" 
{solid line), 2.0" {dashed line), 3.0" {dashed dotted line) and 4.0" 
{dashed three-dotted line) are shown in each plot. 



showing the PA of polarisation. The lowest-level total flux 
contour indicates the level at which 0.8% of uncertainty in 
polarisation is reached. The polarised flux images at J, H 
and Kn are shown in Figure [4] 

At J, we interpret the extended polarisation in the cen- 
tral regions to be dichroic polarisation of starlight through 
aligned dust grains in the host galaxy and the central dust 
lane ( Colina, Sparks, fc Macchetto|199T I . The patchy distri- 
bution of dust in the nuclear regions results in a polarisation 
dependence with aperture size. At H, the most striking fea- 
ture is the polarisation vector pattern (Figures [2| and [3]) and 
the extended structure seen in polarised flux (Figure ET. The 
polarised structure is extended ~3" from NW to SE along 
of a PA ~ 300° . The strong spatial correspondence with the 



structure at 8GHz of a PA of 295° (Morganti, Oosterloo, 
fc Tsvetanov||1998[ ), and the "X-shape" in the [OIII] obser- 
vations at PA of 290° j Colina, Sparks, fc M acchotto| |1991 1 



are shown in Figure [5] The polarised flux images show a 
tentative opening angle of ~75° at the H band, larger than 



the 60° measured through the [OIII] ionised structure by 
[Colina, Sparks, & Macchetto (19911. We note that other 
AGN show a wider opening angle in polarised flux than 
emission line imaging, e.g. NGC1068 has an opening an- 
gle of 60° ±20° from [OIII] imaging (lEvan s et aL]|1991[) and 
80° from NIR polarimetry imaging ( Packham et al.||1997[ ). 
The opening angle in polarised light is larger because of 
scattering can occur from any material in the medium, in- 
stead, the medium needs to be filled by [Oil], in order to 
be ionised ([OIII]), producing a smaller opening angle. We 
interpreted the biconical polarised distribution to be scat- 
tering of light from the central engine to our LOS. At Kn, a 
highly polarised point-like source is observed. 

The polarisation of the central 1.2" of IC5063 at J, 
H and Kn was measured to be 2.0±0.7%, 2.5±0.9% and 
7.8±0.5%, respectively. Figure [T] shows that the degree of 
polarisation decreases as the aperture size increases in all 
filters, whereas the polarised flux remains similar in the 3" 
and 4" apertures at J and H. Finally, the PA of polarisa- 
tion is wavelength-independent (within the error bars) and 
measured to be 3±6° in the three filters. This result is con- 
sistent with the PA of polarisation of 3±2° measured by 
[Hough et al.| ( [T987[ ) in the J, H and K filters, and ~3° using 
optical (0.45 - 0.70 \im) spectropolarimetry by Inglis et al 
(1993'). Note that the PA of polarisation, 3±6°, is perpen- 
dicular to (1) the long-axis, 300°, of the galaxy ([Colina et al 



19911; and (2) the radio-axis, 295° ( Morganti, Oosterloo, fc 
Tsvetanov|l"998| . 



4 ANALYSIS 

The measured polarisation is the diluted polarised light from 
the AGN produced by the diffuse stellar emission. Hence, 
the intrinsic polarisation can be estimated accounting for 
the contributions by the AGN and diffuse stellar emission 
to the total flux in the nucleus of IC5063. We quantify these 
contributions through an examination of the nuclear profiles 
in Section [4. II 
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Figure 2. Total flux contours with polarisation vectors of the central 6" X 6" at J (left) , H {middle) and Kn (right). A vector of 5% of 
polarisation is shown at the top right of each plot. The lowest-level total flux contour represents an uncertainty of 0.8% in the degree of 
polarisation. Contours are plotted in intervals of 10% to the peak of the total flux at each band. North is up and East left. The physical 
scale is 1" = 219 pc. 
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Figure 3. Polarised flux contours with polarisation vectors of the central 6" X 6" at J (left) , H (middle) and Kn (right). A vector of 
5% of polarisation is shown at the top right of each plot. Contours are plotted in intervals of 10% to the peak of the polarised flux at 
each band. North is up and East left. The physical scale is 1" = 219 pc. 
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Figure 4. Polarised flux images of the central 6" X 6" at J (left) , H (middle) and Kn (right) bands. North is up and East is left. The 
physical scale is 1" = 219 pc. 
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Figure 5. Polarised flux at H {grey scale). [OIII] ionisation struc- 
ture from HST by [Sclimitt et aT] (|2003^ {contours). Radio-axis 
with PA~295°is show n ||Morganti et al.||1998^. Du st lane axis 
with PA = 75°, from |Colina, Sparks, fc Macchet"to| | |1991| |. Po- 
larised cones axis (dashed line) with opening angle of 75°it5°are 
from this paper. The physical scale is 1" = 219 pc. 



4.1 The nuclear total flux 

The AGN is embedded within diffuse steUar emission in the 
nuclear bulge. The emission within the nuclear regions of 
IC5063 can be considered to arise from two dominant emis- 
sion components: (1) diffuse stellar emission in the nuclear 
bulge; and (2) emission from the AGN. Photometric cuts 
through the nucleus in the J band show no evidence of a nu- 
clear point-source, and hence we assume that there is negli- 
gible emission from the AGN in our smallest aperture (1.2"; 
~263 pc). Thus, the J band profile is assumed to be rep- 
resentative of the profile of the diffuse stellar emission in 
the nuclear bulge. However, at Kn we detect both a central 
point-source and emission from the diffuse stellar bulge. We 
interpret the central point-source emission to be the AGN. 
To estimate the relative contributions from these two emis- 
sion components, two different methods were followed. In 
the first method, we followed a similar analysis to that of 



Turner et al. (19921 and Packham et al. (19961. We took 



photometric cuts at J and Kn along the dust lane axis, PA 
= 75° ( [Colina, Sparks, fc Macchetto|[l991[ ). These cuts are 
Eissumed to be along a line of constant extinction through 
the nucleus (Av = 0.3 mag; [Colina, Sparks, fc Macchetto| 



(19911). The AGN contribution, modeled as a point-source 



convolved with a Gaussian profile of FWHM equal to that 
of the seeing disc, is termed AGNpgp. To fit the Kn nuclear 
bulge profile along the cut, the emission was modeled as 
the summation of the AGNpgp and the diffuse stellar emis- 
sion from the J band photometric profile. The photometric 
profiles in J, Kn, the modeled AGNpgp and the combined 
J and AGNpgp to fit the Kn-profile are shown in Figure |6j 
showing an acceptable fit to the data. Through this method, 
the best estimate of the contribution of the AGN emission 





Figure 6. Cuts profile for the first method in Section |4.1| are 
shown {lejt panel). Cuts profiles at J {dashed line), Kn {solid line) 
and the scaled AGNpgp {dashed dotted line). The combined J 
profile and AGNpgp matched with the Kn-profile {thick dashed 
three-dotted line). "Flat-top" profile for the second method in Sec- 
tion |4.1| are shown {right panel). Same Kn profile than in left panel 
{solid line) and the subtraction of the AGNpgp-scaled at 45% 
{thick solid line) and 100% {dashed line). 



to the total flux is 40±2% in a 1.2" aperture. This result 
is comparable to that determined for Centaurus A, ~30% 
AGN contribution, when using a similar methodology by 
[Packham eFaL] ( p96l ). 

In the second method, the scaled AGNpgp from the 
first method was used to subtract the AGN emission from 
the Kn image. This procedure produces a "flat-top" pro- 
flle over the center (Figure [6|, which likely leads to a mod- 
est overestimation of the AGN contribution as the underly- 
ing stellar emission will peak at or very close to the AGN 



emission peak (for examples of this technique, see Radom- 



ski et al.[[2002l [2003[ [Ramos Almeida et al.,2009, ^Levenson 



et al. 2009 1. Then, flux of the AGNpgp-subtracted Kn image 
is measured to be 2.25 mjy (13.66 mag.) in a 1.2" aperture. 
Diffuse stellar emission in the nuclear bulge has a contribu- 
tion of 55±3% (45±3% AGN emission) as measured through 
this second method. 

From the methods presented above, the average of the 
AGN emission is estimated to be 43±3%, whereas the contri- 
bution from the diffuse stellar emission in the nuclear bulge 
is estimated to be 57±3%, both estimated in a 1.2" aperture. 
The formal uncertainties are estimated to be 5% from pho- 
tometric measurements. Systemic errors from this method- 
ology could increase the final uncertainties, but are difficult 
to quantify. Previous observations ( [Kulkarni et al.|1998| at 

1.1 ^m, 1.6 [im and 2.2 [im using NICMOS/HST showed an 
unresolved point-source, consistent with an obscured AGN. 
[Kulkarni et al.[ ( [1998[ ) suggested that the infrared emission 
at 2.2 [im is dominated by thermal emission from hot (T = 
720 K) dust in the inner edge of the torus. They suggested 
that the 75% of the flux at 1.6 [im is due to emission other 
than hot dust in the inner edge of the torus, consistent with 
that only diffuse stellar emission contributes to the total flux 
at shorter wavelengths, H band. 

4.2 The nuclear intrinsic polarisation 

The diffuse stellar emission in the nuclear regions of IC5063 
signiflcantly dilutes the observed emission from the AGN, as 
estimated in Section 14.11 We subtracted the diffuse stellar 
emission in the measured degree of polarisation at Kn in 



int 

K„ 



a 1.2" aperture, and an intrinsic polarisation of P 
18.1±1.1% was estimated. Note that the intrinsic polarisa- 
tion calculated here is independent of the polarising mech- 
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anisms of the central engine. The estimated intrinsic po- 
larisation is in good agreement with previous NIR studies 
at J, H and Kn by Hough et al. (19871, where an intrinsic 



polarisation of 17.4±1.3% was calculated, using their own 



photometric data and photometry from Axon et al. ( 1982 1 
in the H and K filters using a 2.25" aperture. Brindle et al. 



1 1990 1 used the same data of Hough et al. ( 1987 \ and with 



slightly different assumptions calculated an intrinsic polar- 
isation of 15.3±1.0% in a 2.25" aperture. They used a po- 
larised power-law component and starlight subtraction. Op- 



tical spectropolarimetric studies by Inglis ct al. ( 1993 1 esti- 
mated an intrinsic polarisation of ~10%, based on narrow 
lines in the range of 0.54 - 0.70 \im. Studies of other AGN 
have shown highly intrinsic polarised nucleus. For example. 



Simpson et al. ( 2002 1 measured a polarisation of 6% in the 



nucleus of NGC1068 and Tadhunter et al. (20001 estimated 



an highly intrinsic polarisation of >28% in the nucleus of 
Cygnus A, where both observations used the 2 [im polarime- 
try mode of NICMOS/HST. 



5 POLARISATION MODEL 

To investigate the polarisation from the torus, the aligned 
dust grains and the role of magnetic fields, we developed a 
polarisation model to take into account the various mech- 
anisms of polarisation in the nuclear regions of IC5063. In 
the sections below, we discuss the possible mechanisms of 
polarisation and then we construct a polarisation model. 



5.1 Possible mechanisms of polarisation 

We consider that the nuclear polarisation could arise 
through three mechanisms in the NIR: (1) synchrotron ra- 
diation, as suggested by Hough et al. ( 1987 1 from a central 



BL Lac type object; (2) dichroic absorption by galactic nu- 
clear dust and the torus; and/or (3) scattering of nuclear 
radiation, as indicated by the ionisation cones. 

5.1.1 Synchrotron emission 

Optical and NIR polarimetric observations of IC5063 by 
[Hough et al.| | |1987p suggested that the large intrinsic de- 
gree of polarisation, 17.4±1.3%, arises from a non-thermal 
nuclear source. The intrinsic polarisation is comparable with 
previous observations of BL Lac objects ( [Angel et al.|197"8{ 
Bailey et al. 1983) and is generally attributed to synchrotron 



emission mechanism. BL Lac objects are observed to show 
a high degree of photometric variability with times-scales of 
weeks and a factor of two in amplitude in the optical UBVI 
filters ( |Impey fc Neugebauer|1988[ ). NIR (J, H and K bands) 
variation, when detect at all, are smaller and longer time- 



scale than in optical (Neugebauer et al. 1989 Hunt et al 



19941. IC5063 does not show any significant fiux- variability 
versus time ( Colina, Sparks, 8z Macchetto|199T I , because we 
are using previous data for calibration, we cannot comment 
further on the issue of variability for IC 5063. However, the 
optical ( |Inglis et al.|1993[ ) and NIR ( [Hough et al.|1987[ this 
paper) polarisation shows a wavelength-dependence on the 
intrinsic polarisation, that is not inconsistent with a syn- 
chrotron source. A further argument against a synchrotron 
origin for the observed NIR polarisation is that the polarised 



radio emission in IC 5063 ( [Morganti e'raL][2007t is located 
1 - 2 SW of the nucleus, with polarisation vectors oriented 
along a very different direction that what we observe in the 
nucleus (Figures[T]and[2|. A decisive test of the synchrotron 
emission interpretation is mm polarisation measurements, 
but as yet no such observations have been published. The 
lack of total fiux variability strongly argues against syn- 
chrotron being a dominant emission mechanism. Hence, we 
do not consider that the polarisation arises dominantly from 
synchrotron mechanism. 



5.1.2 Dichroic absorption 



Young et al. ( 1995 1 showed that the polarisation at K is 



attributed to dichroic absorption of the central engine radi- 
ation passing through the dust within the torus with visual 
extinction 45 mag in NGC1068. They argued that the in- 
crease in the degree of polarisation with wavelength is due to 
dichroic absorption in the NIR. As example, [Simpson et al.[ 
( 2002 1 measured a level of polarisation of 6% in the unre- 



solved nucleus of NGC1068 at 2^m using NICMOS/HST. 
They argued that the polarisation is likely produced by 
dichroic absorption arising from the molecular clouds as- 
sociated with the torus. Hence, the visual extinction of ^ 
45 mag. is consistent with measurements of the degree of 
polarisation of 6%. For detailed discussion of the extinction 
to the nucleus of IC5063 at several wavelengths see Section 
[6] To study how dichroism affects the polarisation in the 
unresolved nucleus of IC5063, this mechanism is included in 
the polarisation model. 



5.1.3 Scattering 

Broad emission lines in polarised flux have been detected in 



the optical spectrum of IC5063 (Inglis et al. 19931, which 
strongly implies the presence of a nuclear dust and/or elec- 
tron scattering screen. Indeed, our observations show, for 
first time in polarised flux in IC5063, the dual-ionisation 
cones at H (Figure |5|, arising from scattering of nuclear ra- 
diation by agents within the ionisation cones. At J, the ion- 
isation cones are not observed as (a) the counter (rear pro- 
jected) cone is obscured by the dusty disc of the host galaxy; 
and (b) the forward (front projected) cone's polarised fiux 
is not clearly detected due to host galaxy contamination 
and/or low signal to noise in the cone area. At Kn, the cones 
are also not observed, as the signal-to-noise of our observa- 
tions is insufficient to map the cones beyond the central 2" 
aperture. To study how scattering affects the polarisation 
in the unresolved nucleus of IC5063, this mechanism is in- 
cluded in the polarising model. 

To understand the polarisation mechanisms responsible 
for the compact nuclear polarisation and the wavelength de- 
pendencies, a polarisation model to simultaneously fit the 
total flux, polarised fiux and the degree and PA of polarisa- 
tion is developed and described in the following sections. 



5.2 Total flux modeling 

To estimate the visual extinction through the torus and con- 
strain the polarisation mechanisms in the nucleus of AGN, 
we developed a polarisation model to fit both the total and 
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Figure 7. Model fit (Section [5]l to ttie total flux in a 1.2" aperture (left panel). The diffuse stellar emission through the nuclear bulge 
{dotted line), AGN through torus {dashed line), ionisation cones {dashed dotted line) and the total {solid line) are shown. Model fit to the 
degree of polarisation {middle panel) and polarised flux {right panel) in a 1.2" aperture are shown. Polarisation produced through the 
dichroic absorption by the dust in the nuclear bulge {dotted line) and within the torus {dashed line), electron scattering {dashed dotted 
line) and total {solid line) are shown. 



polarised flux in a 1.2" (~ 263 pc) aperture of the type 2 
AGN, IC5063. 

The model consists of a point-source emitter, the cen- 
tral engine, partially extinguished by dust within the torus 
and the nuclear bulge. The central engine is described as 
an unpolarised power-law source, (F^ oc i'""), with an as- 
sumed power-law index, a = 1.5, typical of Seyfert 1 nuclei 
and previously used for IC5063 ( [Brindle et al.|1990{ |Colina| 



et al. 1991 Inglis et al. 19931. The power-law is partially 



extinguished by two contributions, one by dust in the torus, 
Av(l), and one by dust in the nuclear bulge, Av(2). We 
also included a contribution from the diffuse stellar emis- 
sion in the nuclear bulge, modeled using typical starlight 
colors of elliptical galaxies, (J-H) = 0.75 mag. and (H-K) — 
0.22 mag. ( Sparks et al.|[l986 l. The diffuse stellar emission 
in the nuclear bulge contribution to the total flux was con- 
strained to be 57% at Kn in a 1.2" aperture (Section 4.1 1. 
The diffuse stellar emission contribution was reddened by 
E(B-V) = 0.26 mag., estimated from the (J-H) and (H-Kn) 
maps of IC5063. Similar values have been used in the liter- 
ature JInglis et al.| ( |1993| l used E(B-V) 
et al.| | |1991[) estimated E(B-V) = 0.1 - 0.4 mag., and |Brindle 
eTaLl ||1990| used E (B-V) = 0.12 mag., calculated by |Bois- 



son fc Durret|(|1986|). A range o f E(B-V) = 0.2-0.6 mag. was 
used by [Bergeron et al. (|1983 K The reddening power- law, 



ta oc A 



by 



Landini et al. 



(19841 was used. 



An additional component added to the total flux to rep- 
resent scattering within the NW and SE ionisation cones was 
also included. Both cones were (initially) assumed to suffer 
the same level of extinction from the nuclear bulge. 

We fit the measured total fluxes in a 1.2" aperture, by 
adjusting the visual extinctions, Av(l) and Av(2), of the to- 
tal flux model. The flt was considered acceptable when the 
deviation from the modeled J, H and Kn total fluxes was 
<1% of the flux value at all wavelengths. Using this proce- 
dure, the fit to the total fiux at J, H and Kn was obtained 
(Figure[7|, giving a visual extinction of Av(l) — 48±2 mag. 



and Av(2) = 6±2 mag. The visual extinction for the host 
galaxy, Av(2), is close to the previous determination, Av = 



7 mag. by Heisler & De Robertis (19991. The total visual 



extinction to the central engine of IC5063 is Av(l) -I- Av(2) 
= 54±4 mag., consistent with the extinction, Av = 64±15 
mag., estimated using optical ([OIII]A5007) and infrared (K 
and L') spectral index and visual extinction maps ( |Simpson 
et al. [1994 1. For the ionisation cones, the SE cone (counter- 
cone) was extinguished with the nuclear bulge visual extinc- 
tion, Av (2) = 6±2 mag, while no extinction for the NW cone 
(forward-cone) was needed. This is entirely consistent with 
the NW cone being in front of the obscuration of the nuclear 
bulge. 



5.3 Polarisation modeling 

Using the same model, the wavelength-dependence of the ob- 
served polarised fiux in the nuclear regions of IC5063 at J, H 
and Kn in a 1.2" aperture was examined. The model assumes 
two separate polarising mechanisms. First, dichroic absorp- 
tion by the interstellar medium of the central engine power- 
law emission. The interstellar polarisation follows a typical 
Serkowsky curve ( [Serkowski, Mathewson, fc Ford|1975[ ): 



P(A) 



exp 



-Kln^ 



/ Amax \ 

\ \ J 



where, Pmax, represents the maximum of polarisation at 
Amax and K = 0.01±0.05 + (1.66±0.99) Amax ( [Whittet et aL 
1992). 

Our polarimetric data and the optical and NIR polari- 



metric data from Hough et al. ( 1987 1 were used to fit the 
Serkowsky curve. The best fit is Pmax = 1.4±0.1%, A,nax ~ 
0.50±0.02 \xm and K = 0.84l°;J*. These values are similar to 
Pmax = 1.49 ±0.08% and Amax = 0.51±0.02 ^m, previously 
estimated bylBrindle et al.l (11990) for IC5063. 



NIR polarimetric studies have shown that the observed 
polarisation in the range of 1.0-2.5 [im is better represented 
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by a power law P oc A with 7 in the range of 1.6-2.0 
(e.g. |Nagata]|1990[ [Martin fc Whittet||1990| ). The degree of 
polarisation arising from dichroic absorption through a dust 
column can be estimated by the knowledge of the extinction, 
Aa, and observed polarisation power-law. Thus, we use the 
expected polarisation by a dust column given by Pa oc Aa 

Second, wavelength-independent electron scattering, P 
= a", is considered to be the polarising mechanism in the 



ionisation cones. Previous studies from UV to NIR (e.g. An 



tonucci fc Miller|1985| [Young et al.|1995| [liiilis et al.|1993 l 

have shown that electron scattering is the dominant polar- 
ising mechanism in the ionisation cones. 

We fit the polarised flux, the degree and PA of polarisa- 
tion using simultaneously the results of the previous section 
with the constraints given by (a) the intrinsic polarisation 
of Pj^* = 18.1±1.1%; (b) the visual extinction of Av(l) = 
48±2 mag., and Av(2) — 6±2 mag.; and (c ) a power-law 
index 7 = 1.9 (Figure [7|. The flt was considered acceptable 
when the difference between the model and observed J, H 
and Kn polarised flux was <1%. Within the 1.2" aperture, 
the polarised flux at J electron scattering (67%) and dichroic 
absorption from dust in the nuclear bulge (33%) are the 
dominant polarising mechanisms, with zero contribution of 
dichroic absorption within the torus. At H, the observed po- 
larised flux is estimated to be 71% dichroic absorption from 
dust in the nuclear bulge, 16% dichroic absorption within 
the torus and 13% electron scattering. At Kn, the observed 
polarised flux produced through dichroic absorption within 
the torus is 89%, 8% electron scattering and 3% dichroic 
absorption from dust in the nuclear bulge. 



Intrinsic polarisation through dichroic 
absorption within the torus 



5.4 



The intrinsic polarisation calculated in Section [4.2| is inde- 
pendent of the polarising mechanisms to the central engine 
of IC5063. To estimate the intrinsic polarisation arising from 
dichroic absorption (P^"^) through the torus, the measured 
polarisation, P^^'' = 7.8±0.5%, at Kn in a 1.2" was cor- 
rected by accounting for (1) the measured degree of polarisa- 
tion through dichroic absorption from dust in an off-nuclear 
region; (2) the estimated polarisation through dichroic ab- 
sorption from dust in the nuclear bulge; (3) the fractional 
contribution of the central engine to the total flux; and (4) 
the fractional contribution of the dichroic absorption within 
the torus. 

The degree of polarisation at Kn produced through 
dichroic absorption from dust in an off-nuclear region of 
IC5063 was measured to be Pg? = 0.2%, with a PA of polar- 
isation ~15°. Note that the difference in PA of polarisation 
at Kn between the nuclear (3±6°) and off-nuclear (~15°), 
represents a lower-limit in the subtraction of the diffuse stel- 
lar emission in the nuclear bulge of IC5063, as polarisation 
is a vector (rather than a scalar) quantity. 

The estimated polarisation produced through dichroic 
absorption from dust in the nuclear bulge is estimated by 
the relationship P(% ) /Av= 1 - 3 % m ag~^ ([Greenberg'l978| 



Whittet et al. 2001 Lazarian 20071. Using our estimated 



mag (Section 5.2 1, transformecj^ to Ak„ ~ 0.7±0.2 mag, 
the dichroic absorption polarisation from dust in the nuclear 
bulge at Kn, is calculated to be in the range of P^'ff = 0.5 - 
2.7%. 

The contributions of the AGN and diffuse stellar emis- 
sion in total flux from Section [4!2[ are 43±3% and 57±3%, re- 
spectively. We define the ratio of the diffuse stellar emission 
and AGN contribution to be AGN™"° = 57±3%/43±3%. 
Through the use of the polarisation model described above, 
the fraction of polarised fiux produced through dichroic ab- 
sorption within the torus is AGNJ^'^ = 89% at Kn. 

The intrinsic polarisation arising from dichroic absorp- 
tion (P^;"^) in the central 1.2" aperture of IC5053 was esti- 
mated as: 



Pdic /T-)Obs 
K„ = l^K„ 



Pk! - Pk5) X (1 + AGN^ 



™"°) X AGN^\^(1) 



The intrinsic polarisation at Kn produced by dichroic 
absorption within the torus is estimated to be P^"^ = 12.5 
± 2.7%. 



6 NUCLEAR EXTINCTION 

The AGN is obscured by the torus and suffers further ex- 
tinction from the host galaxy. Both extinctions were esti- 
mated using the simple NIR polarisation model described 
in Section [5] However, estimates at other wavelengths pro- 
vide further constraints on the extinction to the torus. For 
example. X-ray wavelengths offer perhaps the optimal esti- 
mation of the absorption to the central engine of the AGN, 
whereas the 9.7-(i.m silicate feature in the MIR affords the 
chance to investigate the cooler dust of the torus. In order 
to study the extinction at different wavelengths and its re- 
lation with the torus, the extinction to the nuclear point 
source was estimated through five different techniques, as 
described below. 



(1) X-ray. The visual extinction to the nucleus of 
IC5063 can be calculated using the stan dard Galactic ra - 
tio Av/Nh = 5.23 X 10"^^ mag cm^ (Bohlin et al.[[l978 |. 



[Tazaki et al. (2011 1 calculated an atomic hydrogen column 
density of Nr = 25.0±0.7 x 10^^ cm"^ from Swift/BAT 
data in the energy band of 0.5 - 200 keV. Hence, the visual 
extinction is estimated to be Av = 131±4 mag. A similar 
result of Av = 115±20 mag., was obtained by ^Simpson et alT] 



( [1994 1 using Ginga observations 



(2) NIR polarimetry. Jones ( 1989 1 and Jones, Klebe, & 



Dickey ( 1992 1 examined the relationship between the degree 



of polarisation and the optical depth at K as a function of 
the geometry of the magnetic field. They found that the best 
correlation between the degree of polarisation and the opti- 
cal depth at K for stars extinguished by dust is Pk ~ 2.23 



^3/4 



TjJ (see equation 2 in 



Jones 



19891. This model includes 



with equal contributions from random and constant com- 
ponent of the magnetic field. For the constant component 
of the magnetic field, the relationship between the degree 
of polarisation and the optical depth is Pk = tanh Tp (see 



value of the nuclear bulge visual extinction, Av(2) = 6±2 



equation A7 in [Jones|1989 1. In the case of IC5063, using the 
intrinsic polarisation arising from dichroic absorption within 

^ The conversion of visual extinction, Av, to extinction at Kn, 
Ak„, is Ak„ = 0.112Av l |Jones|1989 l. 
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the torus, P^'' = 12.5 ± 2.7%, and the conversion tk = 0.09 



Av ( Jones|1989 1, the visual extinction is calculated to be in 
the range of 22 - 111 mag. for constant component of the 
magnetic field; and equal contributions from random and 
constant component to the magnetic field, respectively. The 
visual extinction by the torus, Av(l)= 48±2 mag., estimated 
through the polarising model in Section [S] is within the vi- 
sual extinction estimated from these two relationships. 

(3) Silicate absorption. Using the empirical relation be- 
tween the visual extinction and silicate absorption strength, 
Av/rg.T^m = 18.5±1.0, (Whittot 1987) and the sihcate ab- 
sorption, rg.Y^m = 0.32±0.02, measured from N-band spec- 
tra of the nucleus of IC5063 ( [Young et al.|2007| , the visual 
extinction is calculated to be Av = 6±1 mag. This value pre- 
sumably estimates not the extinction to the central engine, 
but the extinction to the outer region of the clumpy torus 
as the observed silicate feature presumably results from a 
mixture of emission and absorption within the torus, this 
should be considered a lower limit to the true extinction. 

(4) Clumpy torus model. The visual extinction to our 
LOS, A^'^'^, is calculated using the clumpy torus model 
i Nenkova et al.||2002| [20081 ) in IC5063. Specifically, the vi- 
sual extinction produced by the torus along our LOS can 



be computed as (see section 4.3.1 in Ramos Almeida et al 



2009): 



A'^"^ = 1.086iVor„e" 



((i-90)V<T^) 



(2) 



where No is the number of clouds along the equatorial di- 
rection; r„ is the optical depth of the individual clouds; i is 
the viewing angle; and a is the torus angular width. 

The most complete fit to IC5063 found in the literature, 
taking into account both photometry and spectroscopic data 
from NIR to MIR and a foreground dust screen geometry, is 
by Alonso-Herrero et al. (20111. We used the probabilistic 



distributions of the free parameters (Alonso-Herrero et al. 
2011 see the blue distributions in their figure A5) in Equa- 
tion [2j to obtain a posterior probabilistic distribution of the 
visual extinction in our LOS (for fu rther details see [Asen 



sio Ramos fc Ramos Almeida||2009| [Ramos Almeida et al. 
2009) . Then, the median and ±lcr from the median were 
estimated to be = 1800^2?!! "ia,g. Note that a flatter 

probabilistic distribution of the free parameters in the visual 
extinction in our LOS, Aj''^, means that a larger error is 
estimated. 

(5) Foreground absorption. The extinction from the dust 
lane in the host galaxy is calculated using the (J-H) vs. 
(H-Kn) color maps of IC5063 and its relation with visual 
extinction ( |Jones|[l989[ ). The visual extinction is measured 
to be Av = 8±2 mag., consistent with Av = 7 mag. obtained 
from dust features in the central 1" - 2" by Heisler & De 
[Robertis] ( |1999[ ); and the value, Av = 6±2 mag., estimated 
using the polarisation model in Section [5] 

A summary of the various extinction values are shown in 
Table |4] The extinction at different wavelengths shows that 
the total fiux obscuration depends strongly on the wave- 
length and/or model used. The estimated visual extinction 
through the dust-to-ratio relation in (1), is due to absorp- 
tion by presumably dust free clouds close to the accretion 
disc, the torus and absorbers in the host galaxy. Instead, 
the estimated visual extinction by the torus model in (4) 
is due to absorption of dust in all clumps along our LOS. 



Table 4. Visual extinction at difTorent wavelengths to the nucleus 
of IC5063 



Method 


Av 




[mag] 


X-Ray 


131±4 


NIR polarisation 


22 - 111 


Polarisation Model 


54±4 


Silicate absorption 


6±1 


Foreground absorption 


8±2 


Clumpy torus model 


1800tf,o 




Homogeneous Toms 



tr-v;; 



Clumpy Toms 



Figure 8. Sketch of the homogeneous {left) and clumpy torus 
{right). For the homogeneous torus, the inner side of the torus 
emits in NIR, the middle regions emits in MIR, and the outer 
side emits in FIR. For the clumpy torus, the NIR emission arises 
from the inner side of the directly illuminated clouds, while the 
MIR and FIR emission is produced from the shadowed side of 
the clumps (see Figure [9]|. Note that clumps located at different 
distances can produce the same emission at a specific wavelength. 
See Section |6| 



Hence, the extinction calculated through these models and 
our polarisation model are produced by different material 
located in and around the torus. Based on these studies, it 
can be interpreted that each wavelength penetrates different 
depths into the obscuring torus. 

The range of obscurations estimated are interpreted in 
terms of the dichroic absorption polarisation. Several NIR 



polarimetric studies (e.g., Jones fc Klebe,1989 Young et al. 
[1995| [Packham et al.[p98[ [Simpson et al.[[2002[ ) of AGN, 
have found that the nuclear polarisation of type 2 AGN typ- 
ically arises through dichroic absorption within the torus to 
our LOS. The implications of this are (1) the total fiux in 
the NIR wavelengths is from the directly illuminated torus 
inner-edge or inner-facing dust clump face (i.e. the surface 
of the dust clump that is directly illuminated by the central 
engine. Figure [s]) or perhaps the central engine itself ( KishTI 
moto et al.|2005[ [2b07, 2008); and (2) the NIR polarisation 
arises as the torus dust grains are aligned, most likely by 
the central engine's magnetic field (see Fig. |9|. 

In this scheme, an individual cloud of the clumpy torus 
can (dependent on the cloud's position and distribution of 
other clouds) absorb radiation from the central engine on 
the inner-facing face. This dust will re-emit radiation at 
NIR wavelengths in all directions. Some of this radiation 
will be self-absorbed by that clump, and some will be emit- 
ted into free space. A small amount of fiux will be emitted 
and penetrate the dust clump at a glancing angle, leading 
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central engine 




single clond 



Figure 9. Sketch of the emission and polarisation of an individ- 
ual directly illuminated cloud in the clumpy torus. (1) The central 
engine emits radiation that is absorbed by the cloud, (2) the ra- 
diation in the outer layers of the cloud is polarised by the passage 
of light through the aligned dust grains, and the NIR polarised 
emission can be in the direction of our LOS. (3) Those rays with 
direction to the nucleus of the cloud, are completely extinguished. 
(4) and (5) the difference of temperature in the cloud produces 
that the warm dust in the middle and back side of the cloud emit 
in MIR and FIR, respectively. See Section [6] 



to obscuration and dichroic polarisation from only a portion 
of the dust clump (see Figure ^ . 



7 MAGNETIC FIELD STRENGTH WITHIN 
THE TORUS 

In this section the magnetic field strength in the NIR re- 
gions of the torus of IC5063 is estimated through three dif- 
ferent methods: (1) paramagnetic alignment; (2) magnetic 
relaxation time; and (3) Chandrasekhar-Fermi method. Al- 
ternatively, and as a possible component of the magnetic 
field strength of the torus, the magnetic field strength at 
Ipc from the SMBH of IC5063 have been estimated (see 
Appendix |A| . 



7.1 Method 1: Polarisation ratio, P/Av(%), vs. 
magnetic field strength 

We attributed the NIR polarisation to dichroic absorption 
from aligned dust grains in the clumps of the torus. The 
alignment can be produced by the rotational dynamics of the 
grain with the environment temperature and/or by the local 
magnetic field ( Davis fc Greenstein||195T l. The orientation 
of dust grains by a magnetic field is called paramagnetic 
alignment, through which grains are become oriented with 
their long axis perpendicular to the magnetic field lines. 

As a first step to considering the magnetic field respon- 
sible for the aligning the dust within the torus, we consider 
the physical conditions and environment of the gas and dust 
within the torus. T he gas temper ature reaches a value of ~ 
10* K in the BLR (|Netzer||l987||. NIR reverberation map- 



ping of several AGN have shown that the outer radius of the 
BLR approximately corresponds to the inner radius of the 
dusty torus ( [Suganuma et al.||2006[ ). [Krolik fc Kriss| ( |2001 1 
suggested that a warm absorber gas in the inner edge of the 
torus can reach temperatures in the range of 10** - 10^ K. 
Recent 3D simulations of the interstellar medium surround- 
ing the central engine showed that atomic gas and ionized 
[CII] trace well the inner regions of the torus, with temper- 
atures in the range of 10'* - 10^ K (iPerez-Beaupuits et al 



Table 5. Physical parameters assumed to the calculation of the 
magnetic field strength in the torus of IC5063 



Description 


Parameter 


Value 


Gas temperature 


Tgas 


10'' K 


Grain temperature 


Tgr 


800 - 1500 K 


Grain size 


a 


10-s - 10-5 cm 


Column density in the cloud 


n 


10"* - lO'"^ cm-3 



limit in the gas temperature to be Tgas ~ W K, in order to 
obtain a lower limit in the estimation of the magnetic field 
strength through methods (1) and (2). The NIR regions are 
located in the di rectly illuminate d clumps in the torus from 
central engin^ [Nenkova et al. (2008.) estimated the dust 
temperatures to be in the range of 800 - 1500 K, for the di- 
rectly illuminated faces of the clumps in the clumpy torus, 
hence we constraint the dust temperature within that range 
in the NIR regions studied in this paper. The temperature 
range is dependent on the size and grain type i.e. graphite 
and/or astronomical-silicates (Nenkova et al. 20081. Based 



on the range of temperatures, the grain size is assumed to 
be in the range of 0.001 - 0.01 (xm. The column densities of 
individual clouds was calculated specifically for IC5063, as- 



suming the following parameters: (1) Alonso-Herrero et al. 
( |2011[ ) obtained a radius of 2.4 pc and a number of clouds 
along the equatorial direction of 14±1 from the clumpy torus 
model; and (2) the gas column densities, as derived from 
NIR molecular hydrogen lines, ranging 1 - 10 x 10^^ cm~'^ 
( Davies et al.||2006[ [Hicks et al.||2009| ). The column density 
for individual clouds in the torus of IC5063 is in the range 
of 10* - 10^ cm . A summary of these physical parameters 
is shown in Tabled 

Models of paramagnetic alignment have been highly 
successfully applied to studies of dust grains in molecular 



clouds (i.e. Lazarian 1995 Gerakines et al. 19951. These 



studies correlate the polarisation, P(%), with dust extinc- 
tion, Av, as a function of the magnetic field strength, B, 



modeled by Vrba, Cyne, & Tapia (19811. The efficiency of 



dust grains, defined as P(%)/Av, is directly proportional 
with the average alignment of the grains, and depends of 
the physical conditions of the environment as well as the 
magnetic field strength, B. An adapted version from equa- 
tion 8 in their paper is presented here to be: 

1/2 



P(%)/Av 



27r 



67x B" 
75aujn \ mnkTgas 



(7-1) 1 



Tgr 



Lgas 



(3) 



where, x", is the imaginary part of the complex electric sus- 
ceptibility, a measure of the attenuation of the wave caused 
by both absorption and scattering; B is the magnetic field 
strength; 7 is the ratio of inertia momentum of the dust 
grains; Tgr, is the grain temperature; a, is the grain size; lj, 
is the orbital frequency of the grains; n, is the column den- 
sity in the cloud; mjj, is the mass of a hydrogen atom; k, is 
the Boltzmann constant; and Tgas, is the gas temperature. 
Purcell ( 1969 1 showed that the lower bound for the most 



2011). Based on these previous studies, we adopt a lower 



3 Note that several authors \Mot et al.|2009||Mor fc Netzer|2012| 
argued that the NIR emission can not be related with dust within 
the torus, but located in inner regions than the inner edge of the 
torus. 
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interstellar grains, the ratio x jio is (see review Aannestad 
[fcPurccll 1973| 



X 



2 X 10"^^ T"^ 



(4) 



The ratio of the moment of inertia of the dust grains, 
7, is defined as: 



7 = 



+ 1 



(5) 



\ is the grain axial ratio. A typical value of hja for 
0.2 ( [Aannestad fc PurceUlfOTS 



where 

interstellar dust grains is 



Kim & Martin 



19951 



Equation ( 3 1 was modeled for optical (V band) wave- 

(|1981|). Further studies 



lengths by Vrba, Cyne, & Tapia 



(Gerakines, Whittet, & Lazarian 



1995 1 showed that it is 



also applicable at K band. Using the physical conditions in 
Table [5] the intrinsic polarisation arising from dichroic ab- 
sorption P^"^ = 12.5±2.7% and the extinction through the 
torus, Av(l) = 48±4 mag, at Kn, Ak„ = 5±2 mag. The 
magnetic field strength was estimated to be in the range of 
12 - 128 mG for the NIR emitting regions of the torus of 
IC5063. 

We can compare our estimated magnetic field strength 
to previous published values. Polarimetric observations with 
the VLA and the GET at 22GHz of the water vapor masers 
in NGC 4258, estimated the value of a toroidal magnetic field 
strength of 90 mG at 0.2 pc (Mo djaz ct al . 2005). Circular 
polarisation observations of NGC4258 estimated an upper- 
limit of the magnetic field strength in the maser features to 
be 300 mG ( [Herrnstein et al.|1998[ ). [Kartje et aI]|T999l ) es- 
timated an lower-limit of the magnetic field strength in the 
maser clouds of AGN to be 20 mG. These studies considered 
outflow winds confined in magnetic field lines generated in 
the central engine. Although for different objects, and at dif- 
ferent spatial locations, our derived range of magnetic field 
strength compares well with these previous studies. 

To estimate the magnetic field strength in the torus 
of IC5063, several assumptions were made. Here, we con- 
sider each of these assumptions in detail. The dust grains 
of the torus clumps almost certainly experience much more 
turbulent and extreme physical conditions than molecular 



clouds (Vrba et al. 1981 Lazarian 1995 Gerakines et al 



1995), which in molecular clouds makes the alignment of 



dust grains less responsive to the magnetic field. 

We assumed a homogeneous magnetic field, where any 
inhomogeneities in the torus are ignored. This assumption 
has some implications for the ratio, P(%)/Av. If a homoge- 
neous magnetic field is responsible for dust grain alignment 
in the torus clumps, then all the dust grains will be aligned 
along the same orientation of the magnetic field line. In this 
case, the ratio P(%)/Av will maximize the alignment effi- 
ciency, and hence the degree of polarisation would decrease 
when inhomogeneities of the magnetic field are present. We 
interpret the magnetic field strength from Method 1, as a 
lower-limit to the magnetic field. 

The method used here shows a strong dependence of 
grain sizes. Based on our interpretation of the extinction and 
polarising modeling, the grain sizes can be constrained. The 
physical conditions in the inner side of the torus (i.e. high 
temperature and direct radiation from the central engine) 
makes impossible any evolutionary grown of the grains in 



such region. Hence, only small grain can survive (0.001-0.1 
(j.m) . Although this physical condition allows us to refine the 
grains sizes ranges, the effect of the grain sizes in the above 
methods are difficult to quantify. 

Thus, our estimate of the magnetic field strength in the 
range of 12 - 128 mG represents a lower-limit for the NIR 
emitting regions of the torus of IC5063. 



7.2 Method 2: Magnetic relaxation time 



Roberge, DeGraff, & Flaherty ( 1993 1 developed a compu- 



tational method to solve the grain alignment problem in 
molecular clouds. They found that the approach followed in 
Method 1 (Section |7.1[ ) is valid only when the ratio of ther- 
mal to magnetic relaxation time is smaller than unity. In 
order to verify this condition, we calculate the lower-limit 
magnetic field strength required to satisfy this condition. 

In the previous section, we assumed that the dust grains 
in the torus clumps are uniquely aligned by the presence of 
a magnetic field. This assumption is valid if the magnetic 
field is strong enough to dominate over both the turbulence 
within the torus clump and the rotational dynamics of the 
clump environment. For rotational dynamic mechanism, the 
dust grains are aligned when the rotational kinetic energy is 
coupled to the grain rotation in equilibrium with the gas 
temperature. The time for the dust grains to be aligned 
by this mechanism is given by the thermal relaxation time 
( |Hildebrand„1988| ), 

1/2 



^thermal 



2ap / TT mjj 



3nmH \8kTgas 



(6) 



where, p, is the grain density. 

In the case of magnetic alignment, the time for the 
grains to be aligned, is given by the magnetic relaxation 
time ( |Hildebrand|1988[ ), 



tmag = 1.6 X 10 



iiaVTgr 

B2 



(7) 



Assuming paramagnetic alignment in the clouds of the 
clumpy torus, then the magnetic field is strong enough to 
align the dust faster than the rotational kinetic energy. In 
other words, the magnetic relaxation time is required to 
be shorter than the thermal relaxation time. i.e. ttij^^^a,! 
> tmag. Using Equations ([g]) and ([7|, a lower-limit of the 
magnetic field strength verifying that condition can be esti- 
mated: 



B^ > 2.4 X 10" a n mnT /^^^Tgas 



Y TT mjj 



1/2 



(8) 



The magnetic field strength is > 2 mG and > 30 mG, for 
the physical conditions shown in Table [5] Note, these values 
are purely theoretical and assume stable physical conditions 
in the clouds of the clumpy torus. Hence, these values are 
considered as a lower-limit magnetic field strength in the 
clumpy torus of IC5063. Since the estimated magnetic field 
strength in Section [TT] are larger than the lower- limit calcu- 
lated here, the condition by [Roberge, DeGraft', fc Flaherty[ 
( 1993]) is satisfied in Method 1. 



Under the condition that the magnetic relaxation time 
is shorter than the thermal relaxation time, the gas and 
dust temperatures are decoupled. This assumption has two 
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implications (1) is only valid in low-density regions in clouds; 
and (2) the ratio P(%) / Av is dependent of the magnetic field 
strength. In our scheme described in Section |6j the detected 
radiation passed through the low-density regions of the torus 
clump, satisfying the condition described above, applied in 
this section. Also, the allowance in the estimation of the 
magnetic field strength in Method 1. 

7.3 Method 3: Chandrasekar-Fermi method 



The Chandrasekhar- Fermi method ( Chandresekhar & Fermi 



1953 1 was also used to estimate the magnetic field strength 
in the torus of IC5063. This method relates the magnetic 
field strength with the dispersion in polarisation angles of 
the constant component of the magnetic field and the projec- 
tion of the mean magnetic field on the plane of the sky (here- 
after termed dispersion of polarisation angles, a) and the ve- 
locity dispersi on of the dust. Here (Equation [9] ), we use the 
equation 8 of Marchwinski, Pavel, & Clemens ( 20121). This 



relationship is an adapted version of equation 7 in |Chan-| 
dresekhar & Fermi (19531 with the factor 0.5 introduced by 



Ostriker, Stone, & Gammie (2001 \ to estimate the magnetic 



field strength in the plane of the sky. 

^4 \ 1/2 



B = 0.5 



[^^G] 



(9) 



where p is the volume mass density in g cm~'^; (j„ is the 
velocity dispersion in cm s~^; and a is the dispersion of 
polarisation angles in radians. 

The volume mass density was calculated using the col- 
umn mass density in Table |5] multiplied by the weight of 
molecular hydrogen. The velocity dispersion in masers ob- 
servations in NGC 3079 is estimated to be 14 km s~^ at a 
distance of ~pc from the central engine ( [Kondratko et ah] 
20051. Based on the magnetohydrodynamical wind model 



for the torus in AGN, Elitzur & Shlosman (20061 assumed 



velocity dispersion in the order of 10 km s~ . Based on these 
previously published results, we use a velocity dispersion of 
10 km s~^. We note that the assumed value of the velocity 
dispersion represent a lower-limit. In order to calculate the 
dispersion of polarisation angles, a, we used the model by 



Jones, Klebe, fc Dickey (1992, hereafter JKD). This model 



relates the degree of polarisation at K with the level of tur- 
bulence in the interstellar medium and the magnetic field. 
In the case of IC5063, using the intrinsic polarisation aris- 



ing from dichroic absorption, Pj,^"^ = 12.5±2.7%, and the 
extinction by the torus, Av(l) = 48±2 mag, we found that 



our data in the JKD model (Figure 10 1 is located between 
(1) the constant component in the magnetic field; and (2) the 
equal contribution of the constant and random components. 
If we assume that the constant component of the magnetic 
field is in the plane of the sky, the dispersion of polarisation 
angles is estimated to be a = 4.5° (0.0785 radians) for our 
data point in Figure [To] In Equation |9] we substituted the 
above numerical values and we estimated a lower-limit of 
the magnetic field strength in the plane of the sky to be 13 
and 41 mG depending of the conditions within the torus of 
IC5063. 

We assumed that the constant component of the mag- 
netic field strength is in the plane of the sky. If the constant 
component of the magnetic field is angled away from the 
plane of the sky, then the magnetic field strength found here 



a. 



JKD Model 










r * * [l^ - 
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^ Gal. ctf. 
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Figure 10. K band degree of polarisation versus optical depth. 
Data from |Jones| l |1989l l {black points) with some of labeled object 
are shown. JKD model assuming constant, = tanhrp, (top 
dashed line) and random, Px t^^^, (bottom dashed line) com- 
ponents; and equal contributions from both constant and random 
^''^ (solid line) are shown. Our data, 
= 48it2 mag. is shown (star). See 



components, = 2.23r^ 

P^''^ = 12.5±2.7% and Av(l) 
Section [73] 



Table 6. Magnetic field strength from three difli'erent methods 



Method 


^min 






(mG) 


(mG) 


1: Polarisation ratio vs. magnetic field strength 


12 


128 


2: Magnetic relaxation time 


2 


30 


3: Chandrasekhar- Fermi method 


13 


41 


4: Magnetic field from central engine at Ipc 


5 


53 



Note: All values represent a lower-limit of the magnetic field. 

will be an underestimate, for example if the magnetic field 
line is pointed to our LOS, we will measure zero polarisation. 

A summary of the estimated magnetic field strength 
through the four methods are shown in Table |6] 



8 CONCLUSIONS 

We presented NIR polarisation at J, H and Kn of the nuclear 
regions of IC5063. The analysis shows a highly polarised 
source, measured to be 7.8±0.5% at Kn, with a wavelength- 
independent PA of polarisation of 3±6° in the three filters. 
For the first time in polarised light of IC5063, the biconical 
ionisation cones are observed, showing a spatial correspon- 
dence with the [OIII] ionisation cones and the radio struc- 
ture at 8GHz, entirely consistent with unified models. 

We developed a polarimetric model to account for the 
various mechanisms of polarisation in the central 1.2" (~263 
pc) aperture of IC5063. To account for the scattering pat- 
tern produced by the biconical structure observed at H, 
an additional polarising component due to electron scat- 
tering was required. The model of the nuclear polarisation 
at Kn is consistent with the polarisation being produced 
through dichroic absorption from aligned dust grains in the 
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clumps of the torus with a visual extinction Av= 48±2 mag. 
by the torus. 

Through the use of various components to the central 
engine of IC5063, we estimated the intrinsic polarisation 
arising from dichroic absorption to be P^"^ = 12.5±2.7% 
at Kn in a 1.2" aperture. Estimates of the extinction to the 
central engine of IC5063 at X-ray, NIR and MIR showed 
a wide variations in the extinction depending on the wave- 
lengths on which the estimated is based on. We interpreted 
as that different wavelengths resulting from different emis- 
sion locations within the torus and hence suffering different 
level of obscuration. In this scheme, an individual cloud of 
the clumpy torus can, depending on the cloud position and 
distribution of other cloud, absorb radiation from the cen- 
tral engine on the inner-facing face. This dust will re-emit 
radiation at NIR wavelengths in all directions. Some of this 
radiation will be self-absorbed by that clump, and some will 
be emitted into free space. A small amount of flux will be 
emitted and penetrate the dust clump at a glancing angle, 
leading to obscuration and dichroic polarisation from only a 
portion of the dust clump. 

We assumed the alignment of dust grains be produced 
by paramagnetic relaxation. Then, the intrinsic polarisation 
and visual extinction ratio, P(%)/Av, is a function of the 
magnetic field strength. We considered the physical con- 
ditions and environments of the gas and dust within the 
torus and we estimated the magnetic field strength in the 
range of 12 - 128 mG for the NIR emitting regions of the 
torus of IC5063. Alternatively, we estimate the magnetic 
field strength in the plane-of-sky using the Chandrasekhar- 
Fermi method. The minimum magnetic field strength in the 
plane-of-sky is estimated to be 13 and 41 mG depending of 
the conditions within the torus of IC5063. 

These studies, to our knowledge, provide the first ap- 
proach investigating the magnetic field of the torus in AGN 
through NIR polarisation. Further NIR polarimetric obser- 
vations of IC5063 and other AGN are required to refine 
and/or modify this approach. The next generation of po- 
larimeters, such as adaptive optics optimized imaging po- 



larimeter in the NIR (1-5 [J.m) MMT-POL ( [Packham fc 
Jones|2010| at th e 6.5m MMT, the MIR polarimeter (7.5-13 
Hm) CanariCam ( [Packham et al.|2005[ ) at the 10.4-m GTG 
will provide a high spatial resolution and polarisation sensi- 
tivity that will allow us to refine and/or modify these stud- 
ies. Also, mm-polarimetric observations with ALMA will al- 
low to refine intrinsic properties of the torus, i.e. dust den- 
sity, grain sizes, temperature, used in the methodology to 
calculate the magnetic fields. 
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APPENDIX A: MAGNETIC FIELD FROM THE 
CENTRAL ENGINE 

The central engine generates a magnetic field. The magnetic 
field at a given distance is proportional to the (a) magnetic 
field of the super massive black hole (SMBH); and (b) dis- 
tance from the SMBIf, follo wing a power-law func tion, given 



byB = B 



BH 



^BH 



(e, 



Silant'ev et al. 



20091. The mag- 



netic field strength at the horizon event is Bbhi a^nd rsHi 
is the radius of the black hole horizon, which is dependent 
of the black hole mass, Mbh- The power-law index, n, is 
taken as 5/4 from the assumed optical thin magnetically 
dominated accretion disc from ftXBH < r < IOOvbh most 



physically significant model of Pariev et al. (20031; and fur- 
ther used by , Silant'ev et al.| ( [2009[ ). In the case of IC5063, 
the black hole mass, 

Mbh = 2.6 x lO'^ M© was estimated 



by Vasudevan et al. (20101, using the mass-luminosity re- 
lation and Two-Micron All-Sky Survey (2M ASS) image at 



K-band. For black hole masses of ~ 10 Mm, Silant'ev et al 



( 2012 1 estimated magnetic field strength at the event hori- 
zon in the range of Bbh = 10* - 10^ G. Using these values 
in the above power-law function, the range of magnetic field 
strength is 5 - 53 mG at Ipc from the horizon event of the 
SMBH of IC5063. 
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